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ieee tT RODUCTION 


A. SOUND PROPAGATION IN SHALLOW WATER CHANNEL 

Experimental investigations of sound propagation in 
shallow water channels have been done by several investi- 
gators. Shallow water propagation is of interest because of 
the applications to coastal defense. These investigations 
are expensive and time consuming. The use of a computer 
model should provide a relatively inexpensive alternative to 
observation. 

One of the techniques uses normal mode theory. The 
normal mode theory, which was introduced and developed by C. 
L. Perkeris [Ref. 1], gave the exact solution in water of 
constant depth. Further development of normal mode theory 
was made by L. Brekovskikh [Ref. 2], who initiated pressure 
as an integral involving Bessel functions and solution of 
the normal mode equation. Another theoretical approach to 
sound propagation in a horizontally stratified ocean of con- 
stant depth 1S given by the method of multiple scattering 
[Ref. 1]. With this method, all the previous theories can 
be simplified by conversion into an asymptotic form which is 
valid when the acoustic wavelength is small compared to the 
distance over which the sound speed varies appreciably. 


These theories agree with the laboratory experiments. 


ll 


For a water channel with a small bottom slope, the sound 
field may be expressed approximately in terms of adiabatic 
normal modes. To facilitate prediction, R.D. Groves, Aqgeem 
Nagl, H. Uberall, and G. L. Zauer [Ref. 3] modeled a wedge- 
shaped isovelocity ocean with a linearly-sloping, perfectly- 
rigid ocean floor uSing adiabatic normal modes. For a pene- 
trable bottom the normal mode description fails when modes 
Oropagating upslope encounter the "critical depth" (H,), 
defined as the depth where the associate mode changes from 
Fully trapped within the water channel to radiating energy 
into the bottom (cut off) [Ref. 4-6]. The parabolic equa- 
tion can be used to explain the mechanism of sound energy 
radiation into the bottom (Ret. 7 palm Such an equation was 
Studied by F. B. Jensen and W. A. Kuperman [Ref. 9], with 
oredictions that satisfactorily agreed with the experimental 
results for small ray angles. With some restrictions, 
normal mode theory 1S applicable for sound propagation in 
the wedge-shaped fluid with a fast bottom. The parabolic 
equation is good for fast and slow bottom, but with the 
restriction that horizontal ray angles must be less than 
ZO 

Another technique introduced to predict the propagation 
of sound in the wedge is the method of images. This method 
was derived from the simplest case; a monofrequency point 
source in a homogeneous ocean with parallel boundaries. The 


total pressure is the sum of an infinite number of spherical 


re 


waves from an infinite set of images. Pie wrest ricei1en Of 
Mercoemetiogmis that it sdoes not. generalize to the case of 
inhomogeneous media or non-planar boundaries. In this work, 


mrrcemethod wili be studied. 


pee FHE MNEFHOD OF IMAGES 

hilo voeNCODpens, cancers, [Oannou, and Kawamura [Pef. 
10], predicted the pressure amplitude and phase of the sound 
field along the bottom of a wedqe-shaped fluid layer of 
Meicity pi, and speed of sound cj), overlying a fast fluid 
Berto of density o>, and speed of sound c> > c, by applying 
the method of images in a computer program implementation. 
In 1984, Baek [{Ref. 11], and LeSesne [Ref. 12], implemented 
further improvements. Baek's computer program, WEDGE, and 
LeSesne's computer prodqram XSLOPE were validated for several 
cases. WEDGE was developed for two-dimensional upslope pro- 
NMaqation (the source and received are in the same vertical 
nlane pernendicular to the shore line, and the receiver is 
Meescer tO the apex than the source (Fiqure l.la)) and down- 
Slope pronaaqation (the source is closer to the apex than the 
Beeeivyer (Figure 1l.l1b)). XSLOPE was develoned for upslope, 
downslope, or cross-slope propaqation (the source, receiver, 
and apex, are not necessarily in the same plane nerpendicu- 
MmmeetOo the shore line (Fiqure 1l.lc)). [Reborn MpEoag rans, 
Raek and LaSesne assume that the fluid in the wedqe and 


mard = §8iNn hic Ot tom have eonsitea tle densities, teats 


ee 


the speed of sound is constant, and that the  intertacemaem 
tween the fluids and the surface is smooth. 

In both WEDGE and XSLOPE, all distances are "scaleqmeme 
UNLES "OL tee Gdunpedlst ances A dump distance X, as stated 
in Reference 10, is the distance from the apex measured 
along the interface at which the lowest mode attains cutoff. 


If the wedge angle is 8 (Figure 1.lc), then 


cms m2 


kK} Sin, tang (1.1) 
Gens CbeGOs (ceo. (loa 


where k,; is the wave number in the wedge and 8, is the Cries 


ical grazing angle for reflection of sound from the bottom 


FOres << J 
X = H tan 8 (lon 


This scaling distance negates the necessity of specifying 


frequency. 


ce COMPUTER PROGRAM DSLOW 

At the start of the work reported in this thesis, a com- 
puter program was obtained [Ref. 13], which is an extension 
of the WEDGE and XSLOPE for downslope configuration with a 
slow bottom. The computer model, DSLOW, developed to run on 
a desktop computer (Wang 2000), uses the method of images to 


predict the pressure amplitude and phase anywhere 
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within the wedge fluid overlying a slow bottom in a cross- 
EitepercOntiaguration. A geometrical picture of this configu- 
motion 15 Shown in Figure l.lc. 

Mathematically, the model used in WEDGE and XSLOPE is 
meplicable in any condition. But consideration must be 
given for making 1t work for a slow bottom. In the case of 
a fast bottom, the dump distance has a physical meaning. 
The dump distance is expressed as a function of the critical 
anqle. The critical angle is equal to arccos (cj/c9). In 
mie case Of slow bottom, cj/C> 1S greater than 1, thus 
Pieecos (c]/C>5) 1S invalid; therefore, so is the dump 
distance. To facilitate the scaling factor, a "characteris- 
tic distance" or "scaling distance" is introduced. We need 
the scaling distance because, with this distance, our model 
will be independent of frequency as in the fast bottom case. 
There is also the hope that the use of a scaling distance 
will allow systematic observation of the pressure field. 
Mims, Scale distance X, 1s the distance measured alonq the 
interface from the apex to the point where the lowest mode 
meulad attain cutoff if the fluids in the wedge and in the 
bottom were to be interchanged. The characteristic distance 


is defined by the following equation: 


_ 1/2 
Sona k> siné, tanB (1.4) 


LS 


where 


80 


= areees (C77 cam and K5 = w/C > is the wave number 


lueetne DSoOtLLom. 


ie 


The following terms will he used throughout (see Figure 

B = weccre anqle 

R; = distance of the source from the apex in) Sua 
Xo 

R> = distance of the receiver from the apex in (UinaieSaaam 
x 

iv = angle of elevation of the source above the bottom 

§ = ancle of elevation of the receiver above the 
elena wie. 

Yo. = Gistance bere cniare projection of the source and 
receiver on the shore line veseabecdmov.«. 

p]/e7 1s the ratio between the density of the — Saaismamy 
the wedae (p,) and the density of the fluid ingen. 
bottom (99) 

c}/Cz 1s the ratio between the Speed of sound)” (ae 


wedge (c,) and the speed of sound in the 5H@tiiam 
(eo A fast bottom occurs when C5 > ¢);7 a Siam 
DOttom occurs when Gc) ae) 


The nurpose of this research is the following: 


1. To transfer, test, and evaluate DSLOW program On ieee 


Lew 33 00- 


2. To obtain numerical and qraphical out puta eee 


number of cases; and 


3. To attempt to develop plausible explanations for any 


Slanificant features observed. 
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fa Upessope Configuration 
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b. downslope configuration 





6. eross-slope configuration 


Figure 1.1 Geometry of the wedge 
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ll. THEORY 


A. GENERAL VIEW OF A WEDGE PRESSURE DISTRIBUTION EN > THE 
DOWNSLOPE CONFIGURATION 


A general picture of the sound energy propagation within 
the wedge in downslope direction is given in Figure 2.1. If 
a sound source is placed at point S, ray 1 will reach the 
Surface at point P with an incident angle a with respect to 
the normal to the surface at this point. This ray 1s re- 
flected by the surface at the same angle but with the phase 
P80? “dil ererent. (On the surface, sound pressure 1S zero 
everywhere.) The reflected ray reaches the bottom with an 
Incident angle 8 +a. At great enough distance, ray l 
never reaches the bottom again. This ray does not contri- 
bute to a sound pressure field at the bottom. The pressure 
at the bottom should be very small according to the ray 
theory argument. 

Using these ray-tracing methods, an estimated profile of 
the pressure amplitude versus the receiver depth can be 
made. When the source and the receiver are placed near the 
apex, the pressure amplitude is zero at the surface, a maxi- 
mum somewhere within the wedge, and greater than zero at the 
BOEEO: In the case where the source is at a far distance, 
the pressure amplitude is equal to zero at the surface, a 


maximum somewhere within the wedge, and zero at the bottom. 
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Payeeeracuna will only give a rough approximation, not an 
exact solution, but ray tracing may be used as a guide. The 
method of images calculates the exact pressure amplitude at 
each point within the wedge subject only to the assumption 
inherent in uSing the plane-wave Rayleigh reflection 
coefficients. 

B. SOUND PRESSURE AT A POINT IN THE WEDGE DOWNSLOPE 

PREDICTED BY THE METHOD OF IMAGES 

Let the source be a scaled distance R, from the apex and 
at an angle of y meaSured from the bottom of the wedge. Let 
the receiver be a scaled distance Ry from the apex and at an 
angle § measured from the bottom. 

Using Figure 2.2, let the upper half family of images be 
n = 1,2,3,4,... and the lower half family be n' = 
ie, o,4..-. Calculating the field resulting from source and 
images proceeds along the lines developed in [Ref. 14]. it 
meee ene angle stormed at the apex between the ny, image of 


the source and the receiver, then 


pea eo aay 
Ce eB On ageny 
oe) See 
eee et ee = Oraceenty 
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Ray tracing in the wedge downslope 


Figure 2.1 
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Or more generally 
few net Bes 16 > 7 Lone Todd 
esas 0 te mo vaeeile ey. CN 


Which can be reduced to: 


rear 2 l= ( =) eel) ly - 6 (20) 


or 


n 1 


+ 
bn 2 iNT | 5 }B + (-1)"y - § (22) 
where INT[{ ] denotes the largest integer which is equal to, 
or smaller than the argument. Using the same method for the 


member n' of the lower family of images we obtain: 


bn! {n+(1/2)[1-(-1)"}}g8 + (-1)"y + 6 (2 moe) 


n 1 


bad = 2 aN s—— 18 + (-1)"y + 6 (2.4) 


Using the geometry of Figures 2.2 and 2.4, the distance 


eA th 


between the n and on images to the receiver is 
resnectively 
ee ¢ R}- + R>- - 2R )R5cC0S6, (Ze) 
and 
Age Seo) oe ae eee 
rn = YR]4 + Ro* - 2R]R7C0S$9,, (236) 


Zu 


ee Sa | 


je. 


QO image = source point 


Figure 2.2 Geometry of a wedge by the method of images 


Lee 


iemeaiGmes sO and Goae for the nth ang ntth images 


resnectively are 
SIn oy 


8no = arc tan [ Ro/P] - cos 6, ] (2 


and 
SIN¢y,! 


Onto = Aas tan | RSVR] > GOS One | (2 08 


ae ec otro che angles of incidence for the 


m=h bounces from the bottom for the n and n' image 
Mmrmectively:; m= 1,2,3..6.«. (The 0&2 bounce is the last 
one before reaching the receiver.) The geometry of Figures 


feeeand 2.3 aive 6,, as follows: 


65) = 899 - 28 - 8 
639 = 939 - 48 - «6 
one | 
Seema 0. as 


The qeneral expression is 


ee 2mp - 6 


Using the same method 


= @.,:. - 2mp + 6 


On'm Oete: 


fei (ae 


mhe Maximum numoer of bottom bounces of the n ane nn 


image is 
eee es IN lo 726) = INT [¢,1/28] 


Ze 


apex 





Figure 2.3 Geometric development of reflection angles 
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The maximum number of images is 


Bees = ON = INT {17 8) 


The reflection coefficients for the n&® ang ni th images 


for a plane wave are: 


po] 


- 
_ ~2°2 va aT 
X39) 9 ————————— or 
GC + Ynm 
P2= 2 
and 
pP1o] _ Ye 
p2&2 (2. UG) 
Raper.) 1 
+ VUn'm 
P2092 
where 
Perel Ce iecoss mo. 2. 
Ynm S neo (2.11) 
and 
mm, = Ce COS a 0 ie (22) 
¥n m sin On 'm 


The contribution from the upper family of images is 
M 
I (253) 





N 
Pu = ) : exp(-jkr,)(-1)INTIn+1)/2] 
=] ‘n 


Rnm 
n =(0 


a 
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Figure 2.4 Geometry of symmetric images 
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and for lower family of images is 


N 
Pu = } = exp(-jkr,1)(-1L)INTEnt1)/2] eae (2.14) 
n=] —9 m= 0 


The total complex pressure is 


P(x) = Pu + Pl Vo) 


oy), 


III. DSLOW PROGRAM IMPLEMENTATION 


A. PROGRAMS FEATURES 

Since the mainframe graphics computer was available, the 
DISSPLA graphical program was used. The only programming 
language compatible with DISSPLA is FORTRAN. The numerical 
and graphical output is provided by this program. To give 
the pressure amplitude versus received angle graphs, two- 
dimensional plotting 1s used. 

The program DSLOW is run by placing the point source 
anywhere in the wedge and then placing the receiver at a 
distance downslope from the source. The receiver position 
was varied from zero degrees at the bottom to 8 at the 
Surface. High resolution plotting was achieved by dividing 
the y-axis (received angle) into two regions. The filmer 
region covers the receiver angles from zero to 1/5 of the 
wedge angle. In this region Aé is equal to 8/100. The 
second region covers the remaining wedge angle with Aé equal 
Co 8, LO. This method provides 29 predictions of the pres- 
Sure amplitude. Another method of plotting carried out was 
in the region of 6 > 8/2, Ad = 8/10, and in the region Jot 
< B/2, A&= B/100. This method provides 54 points to be 


plotted. 
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Br NORMALIZATION 

The main goal of this research was to investigate the 
profile of the pressure amplitude as a function of a number 
of variables. An example of the numerical values of the 
pressure amplitude, the normalized pressure amplitude, and 
the phase at each receiver position is displayed in Appendix 
ec; The sound pressure becomes smaller as the receiver is 
moved away from the source. It the pressure amplitude were 
plotted directly, it would be difficult to compare the 
Curves at near distances to the curves at far distances 


Since at the near distances the pressure amplitude 1s much 


greater than the pressure amplitude in far distance. Thus, 
a normalized nressure amplitude is needed. The normalized 
Mressure 1S obtained as follows: (see Figure 3.1) 


We know that the sound pressure at the surface 1S zero 
and that the sound pressure is a small number’ greater than 
zero at a point near the surface. The first non-zero value 


mmmperessure P) 1S at the receiver angle, 6, = 98/10. We use 
this first calculated non-zero pressure amplitude as the 


normalization unit. The normalized pressure is 
PN = P(6)/(P}) (2) 


where P(6é) 1S the pressure at any point within the wedge. 
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receiver angle 


. curve of the 
NG normalized pressure PN 
q—-————_————. VS. receiver angle 


- 


curve of the 
true pressure P (64 ) 
vs. receiver angle 


oi pressure amplitude ——», 


Figure 3.1 Pressure amplitude normalization 
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oe PROCEDURE 

Figures 3.2 through 3.7 represent the results when the 
receiver distance and source angle are fixed and the source 
distance and receiver angle are varied. Figures 3.8 through 
3.13 renresent the results when the source distance and 
angle are fixed and the received distance and angle are 
varied. These cases will be the foundation of onas 
subsequent discussions. 

The solid lines indicate the fitted curve and the dots 
indicate some values of the normalized pressure amplitude. 


In DSLOW, the dot appears at each third datum. 


Dr Peoscrat IMPROVEMENT 

DSLOW was designed to provide three-dimensional graphs. 
For example, the x-axis represents’7 the scaled source 
distance, the y-axls represents the scaled received 
distance, and z-axis represents the normalized pressure 
amplitude. Lo samo) £y the presentation, only two- 
dimensional graphs were presented with the x-axis the nor- 
malized pressure amplitude and the y-axis the receiver angle 
6. All curves are presented with the data fitted with a 
cubic spline. 

The DSLOW program was executed to obtain numerical 
results of the phase angle, the pressure amplitude, and the 
normalized pressure amplitude at each receiver position. 
The iminSte run used double MEeC1SiLon ine eae Clibacy-. 
Difficulties were encountered when the DISSPLA subprogram 


Se 


was attached for making the graphical output. When double 
precision and DISSPLA were not successful, the single pre- 
cision was used, resulting in round-off error. (See Figure 


3.8 at Rp = 10.0.) 
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Figure 3.6 
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amplitude with Ry fixed, R, varied 


Figure 3.9 
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amplitude with R»y fixed, R, varied 
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ive DISCUSS LONS 


eee GRAPHICAL OUTPUT 

The graphs of normalized pressure amplitude as a func- 
tion of receiver angle 6 were investigated for various 
Powmee distances Ry and receiver distances R>5 (Fig. 3.2 to 
3.13) while the other parameters are held contstant. Ona 
given wedge angle 8 and sufficiently small source distance 
(Fig. 3.8 and 3.12), at all receiver distances, the pressure 
increases uniformly towards the bottom. For greater source 
distances, (Fig. 3.9 ~-3.11) the pressure attains a maximum 
meenin the wedge for all receiver distance. 

As the receiver distance 1S increased (Fig. 3.9), a 
pressure minimum develops between the maximum and the 
PoOEtom. An important property of the curves of pressure 
versus receiver angle when there 1S a maximum and minimum is 
that, at a specific receiver distance, the pressure above 


the minimum can he extrapolated to zero pressure on the 


Bottom. foccmrigu, 3.10 wither, = 32.)) This receiver dis- 
Mmamece 1S Called the “transition point." So far, we do not 
mew the properties of the transition point. We use the 


transition point for indicating the behavior of the curves 
when the parameter involved 1S varied. The transition point 
appeared twice in some cases, but in the following 


discussions the first transition point is the only point we 
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will be concerned with. (See Fig. 4.1 for transitions Cores 


respond to Ry = 4.6 and 6.4.) 


B. GRAPHS CLASSIFICATION 

The development of curves with the source distance (R})) 
and the receiver distance (R32) as variables was observed. 
As R,} or Ro are varied the curve changes from a line@ap@euma 
to a curve with an observable minimum (Fig. 3.9, Rza = 5.9) 


and finally to a curve without a minimum (Fig. 3.9, 3a 


a0, Three different types of curves resulted from the 
series of two-dimensional plotting. They are described 
below: 


1. Type l Curves 
Type 1 curves (Fig. 4.2) are those where the sound 
pressure 1S equal to zero at the surface and maximum at the 
bottom and 1s almost linearly dependent on depth. These 
Curves are most pronounced when the source distance is much 
smaller than the characteristic distance. The closer’ the 
source 1s to the characteristic distance, the more nonlinear 
the curves (Figs. 343° and 3405. 
2. Type 2 Curves 
Type 2 curves (Fig. 4.3) are those where the sound 
pressure iS zero at the surface, maximum somewhere between 
the surface and the bottom with no minimum. These types of 
Curves are generated when the source is placed at a point 


much greater than the characteristic distance. Type 2 
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Figure 4.1 The plots where there are two transition points, 
R9=4.6 is the first and Oo acwee Ss tne, second 
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Curves indicate that the sound energy in the wedge is well 
collimated and that reflection is negligable. 
347 1Y Deno Curves 

Type 3 curves (Fig. 4.4) are those that have a min- 
imum pressure. These curves occur when the source is a dis- 
tance slightly greater than, or less than, the charac- 
ter rstl1c distance. Tables 1, 2, and 3 of Appendix OD show 
the receiver positions at the first transition poumeee 
Three different values of 8, two different values of pj/P9, 
and two different values of c)/c j were used in making these 
tables. The transition point did not occur when 8 = ]5e8 
pee Gon») On 07 ea coe Oe An explanation can be offered 
using the fact that for these particular sound-speed and 
density ratios an angle of intromission exists [Ref. 15]. 
Since the angle of intromission 1s the grazing angle at 
which the sound energy 1s completely transmitted into the 
slow bottom, it 1s plausible that no transition “pou 


GeCursS. 


Ce JTRANST TION. 2OLN a 

By varying the wedge angle 8 in small increments Ag = 
0.5° starting with g = 5°, and ending at 8 =" 7 eee 
found that transition occurs for source distances within the 
range £rom 170 Prowl. o 

For Ry < 1.0, no transition point was observed; )eem- 


curves are the Type l. For 1.0 < Ry < 1.5, the evolutionges 
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curves as the receiver distance varied can be explained as 
follows: first, the receiver 1S Dlaced near the source and 
gradually it is shifted further from the source. The mini- 
mum in the pressure decreases reaching the point where tne 
Curves extrapolate to zero (the first transition point). 
Further detailed observations were made on this particular 
Facet by varying the source distance and the receiver 
distances. The results of these observatini1s are tabulated 
arc craphed in Appendices D and E. When the recelver 1s 
moved away from the source, the minimum will reach a minimum 
pressure then the pressure increases until it reaches the 
point where the curves again can be extrapolated to zero, 
this is the second transition point (See Fig. 4.1). 

BOmery 2 l.5, there will™ge no transition point. The 
curves are the Type 2. 

The transition point as a function of source angle can 
be observed using the tables in Appendix D. In most cases 
the greater the source angle y, the closer the transition 
point is to the apex. Graphs of transition point as a func 
Mmeoren, (Appendix EB) indicate that the smaller 8 the more 
regular the curves. This is easy to understand because the 
smaller 8, the more accurate the observation of transition 


Bernt; the greater 8 the less accurate the data. 
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Type 2 curves, indicating a pressure 


amplitude nearly linear with depth 


Figure 4.2 
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Figure 4.4 Type 3 curves, indicates the presence of 
reflection and refraction near the bottom 
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Dy PARAMETER VARIATIONS 

Variation of parameters was done by changing one 
parameter of interest while all others were held constant, 
for fixed source and receiver distances, and plotting the 
receiver angle versus normalized pressure amplitude. 

The parameters 8, p}/P9r and c)/cy, were held constant 
and the pressure amplitude was plotted for various 56, Rj, 
and R>- Tiemmarsmane Set at 6/4, 6/2, and 36/4. Variations 
momeie Shore distance (Y,) can be made because the program 
is available, but to simplify the investigation, Y,y was set 
emai to zero for all plots (Fig 4.2 is included as an 
example for Yo # 0). 

if) Velabatioms of 6 

Initially, the values of g investigated were: O% 
me and 15°. The major effect created by altering the 
value of g is that, for the same values of Rj, 1/2, and 
peo, ene Smaller 8, the shorter the transition point (see 


pei 3.9? and 3.10, Tables 1, 2, and 3 of Appendix 0D). 


Poe ariations@or y¥ 
The Variations of y from y = 6/4 (the source is 
placed near the bottom) to y= 38/4 (the source is placed 


near the surface) are presented in Tables 1, 2, and 3 indi- 
cated that the greater y the shorter the transition point. 
Mmmeoemotc always true, for instance in Table l at 6 = 6°, Rj 


=el.o0, the greater y the longer transition point, for the 


a3 


rises. (See Appendix D.) 


3. Variations of ¢)/¢c> anda a. 
Variations of the acoustical parameters c)/c 9 and 


p1}/e7 were done, but did not give a significant variation of 


the sound pressure profile. 
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Ve CONCLUSIONS eS RECOMMENDATIONS 


fee CONCLUSIONS 

The program DSLOW gives good plots representative of the 
sound energy distribution within the wedge. The sound energy 
can be well collimated by the wedge. This phenomenon 1S 
eeeongly aftected by the source position. AG “a “source 
NMOsition close enough to the apex, sound energy is distri- 
buted linearly with respect to the depth. As the source 
moved away from the apex, the distribution of sound energy 
becomes more complex. Sometimes a minimum is7~ found; this 
minimum may be caused by the presence of sound energy 
reflected by the bottom. 

The source position plays a major role in forming the 
pressure distribution profile. The pressure distribution is 
also very sensitive to the parameter variation at small 
source distances, but it becomes insensitive at large source 
distance. The characteristic distance must have physical 
meanings rather than just an arbitrary number, because when 
the source distance in proximity to the characteristic 
distance, the model is most sensitive. 

The model is restricted when the single precision mode 
generates round-off error and rough curves which do not 


allow for accurate analysis. 
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Bae 


RECOMMENDATIONS 


ee 


Single precision produces good results, but fang 
in some cases. Double precision would improve the 
program, but increase the execution time. This emer 
be done by running the program in double precision 
and accumulating the result in single precision 
before plotting the data by DISSPLA. 

[It is sugqaested that the program be run using more 
realistic narameters and observing the effects on 
the characteristic distance and transition point. 
Further study validating DSLOW in comparison with 


exnerimental results 1s suqgested. 
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APPENDIX A 
DSLOW ALGORITHM 


The pressure amplitude calculation 
Rp LNT (1S0'p) 


r 
Pax 


 2tanB tan[arccos(c>Cc))] 
AL=a'k,=0.0001 (constant) 


me 2 2 2 
Pmt g TRY * Rp 
Rg=V D5 - R3 cos((N-1)B+y - 4] 
Rg= JD» - R3 cos((N-1)B+ y+ 4] 
$5 =(-1)ESTOSy 2) 
W,=2AL(cy/c5)* 

] 2 


SI=|{Ry sinf((N-1)B + y] - 2ENT((N-1)/2)B 


+ Ry sinf2INT(N-1)/2)B - SPV Rg 
Cl= J(1-SI7) 
T= SI’D, 


ben (o> C]-C>) 
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(eoneAya) 


(eqn A.2) 


(eqn A.3) 


(eqn A.4) 


(eqn A.5) 


(eqn A.6) 


(eqn A.7) 


(eqn A.8) 


(eqn A.9) 


(eqn A.10) 


(eqn A.11) 


(eqn A.12) 


(eqn A.13) 


Y=\ae: Wc 


Ca a 
Yy = wy YG 


Y5=V(¥-W,),2 


=. ele 
TNS 
2 l 


aa 


a TS 
‘~ Ci ea 


(T?-¥5) = Na 
= 
a eN5) so 


2¥)T 
26 ne 
(TVs) +) 
N 
P= © iyi 124. cos(Ro_ Kk X)+ Ze sine ene e em 
fiat 
.) 
Pa= ¥ (-1 STO 2¢-Z. sin(Rg K, x) + Ze eosin a Ree 
n=1 
N 


err 
P3 =e el)! CN 2)(2. cos(Re_ Ky X)it Zeist Rema aa 
(Via 
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(eqn A.14) 


(eqn Avia 


(eqn A.16) 


(eqn A.17) 


(eqn A.18) 


(eqn A.19) 


(eqn A.20) 


(eqn A.21) 


(eqn A.22) 


(eqn A.23) 


(eqn A.24) 


< 


Be N, 
Pye Y (-1 TN 2)1.Z, sin(Rg Ky X)+ Zp cos(Rg pK pX}/Rop 


n= 1 
Gage Py 


2 aes 
Brot a RIV Ps” + Pe 


a0 


(eqn A.25) 


(eqn A.26) 


(eqn A.27) 


(eqn A.28) 


APPENDIX B 
DSLOW PROGRAM 


FI I IK I IKK KKK IKK KKH KR KKK RR KK AKER KEK K KKK KKK KKK KAKA REKKKKKKKKKK 


ae THIS PROGRAM CALLED DSLOW IS_ CALCULATING THE SOUND PRESSURE ne 
a WITHIN THE WEDGE OVERLYING SLOW BOTTOM FLUID AND DOWN-SLOPE on 


REC 


DO eae eee aaa eee BO siete ic cik skied eR A 


INTEG R A,I,11 MN Si s2.Niedere 
2 B. CC, 22 gp Pt. b2,G, Piociece QL RL R2 T, 


i oF PS, TQQ4 6422 Hi i gto) 8 1365,5(35) 5730), E1409 3 ba bn2 46) 


RARAAKRRAR AEDS =4,9 9090 


INPUT PARAMETE 


Be RS rede RI TOI II OA I AOR ea ea rea else ok se ech fe see eS 


a 


KAAEKRKEKEKKERRKERKERKRKEKRKEKKREKRRKKKKKRKKEKKKKEKKEKEKKKKKKKAKK 
x* 


Dee xx 

Ee B = WEDGE ANGLE ( oe 

C G = SOURCE ANGLE é 

G D = Sea ANGL DEG) 

C N1l= OF IMAGE POINTS 

C Me = SOURCE DISTANCE {i CHARACTERISTIC DISTANCES 2 

Cc R2 = RECEIVER DISTANCE: (IN CHARACTERISTIC DISTANCES) 
Cc YO = APEX DISTANCE (IN CHARACTERISTIC DISTANCES 

C Dl = RHO 1/RHO 2 

CG CC = C-1/C-2 

g AL = ao bE4 ITIO 

C ern tan nn ee Qe RECEIVER POSITIO ND se dedcdk dead dee dee de ek dea ek ak A Seok eek A 
K* 

SIC eee ee eee ONS ekkkkkddha hank kkk K RRAK AAR ARR 
G YO = 0.00 

- RZ = 9.00 

- DO 111 P=1,6 

C B = 10.0 

- B = 10.0 

C G = 3*B/4 

c G = 2*B/4 

> G = 1*B/4 

- D= 0.0 

C Rl = 0.95 

- Ri i=* Leto 

G Rl = 5.99 

- Rl = 0.7? 

c RZ = Ok 

E R2 = 1:2 

C R2 = 1.50 

C oS = 2.00 

g RZ = 5.00 

- TU = 0.00 

- D1 = 0.80 

- CC = 4-48 

- AL = 0:0001 

Lede ek ht REAR URKHHARK IK RK HHH K KAI KRM 
x* x* 
FO ee Ne Coe eee eR ee ee oe ROSS aka ane eK 
- K = Q 

- Nl = IN ( 80./B) 

- T6 =_18 PI 

- B = B/T 

- G_ = G/T 

- C2 = CC**?2 

- bg = (¥Q*¥Q)4(R1*R1)+(R2*R2) 

g T4 = Pill * ae 

- i An B 

- i 3 = AC il 

- ly = TAN 

- TQOO3 = x 002" 

- T4*= PI/TOO 
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e CALL MESSAG('SHO,DIST: Ss 
€ CALL MESSA + Bt 068", 100, 'AB ut’: (ABUT » 
C CALL RESET('ALL' 

C CALL PARA3 

C CALE POLY3 

€ ERLE CURVE (PN, DZ, 29,0) 

C111 CONTINUE 

C CALL ENDPL(O) 

G CALL DONEP 

C SDOF 

C END 
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APPENDIX D 
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Pen Ra ae ae a 
y= | 324 | 376 | 47s] 73 | 8 


B=10°, K,X=19.44 


| Rye] RH tao | RpHisd | Ry= 1] Kya LO 

py=Ba | 98 | uo [| io | 400 | 600 | 

py-B2 | 76 | 36 | to | 190 {400 
= | : . QI = 











tes 


a 


APPENDIX E 
GRAPHS OF R, VERSUS R, AT THE FIRST TRANSITION POINT 
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R1! VERSUS R2 AT TRANSITION POINT 
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Figure E.1 Ry vs Ry at the first trans.points, for B= 6°, P ;/Pz = 0.80, cp /cy = 1.10. 
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0.80, Cy/C> = 1.10. 


Figure E.2 Ry vs R> at the first trans. points, for B= 10°, Py /P> 
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Cy/C4= 1.10. 


0.80, 


Figure E.3 Ry vs R» at the first trans. points, for B= 15°, 91 /P> 
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Figure E.4 Ry, vs R» at the first trans. points, for B= 6°, P 1 /Pz = 0.80, cy/c4= 20. 
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C1 /C> = 1.20. 


Figure E.5 Ry vs R> at the first trans. points, for B= 10°, py/p>= 0.80, 
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Figure E.6 Ry, vs R» at the first trans. points, for f= 152. Pp )/P> = 9.80, c)/cy = 1.20. 


2 


2° 1=29/19'9'=ZOHY/10HY'(9490)9=4 
yi 
eT Q's Lt 9°] oT rt 


VLAG=VNAV) _o. 
6/VLAEG=VAWNVD o 
QNADgT 


ENIOd NOILISNVUL LV cod SASHA Ta 


Cl 





Figure E.7 Ry, vs Ro at the first trans. points, for B= 6°, P )/Pz= 0.90, cy/c2= 1.10. 
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Figure E.8 Ry, vs R> at the first trans. points, for B= 10°, P )/P> = 90.90, c)/cy= 1.10. 
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